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Abstract. We investigate the kernels of the transformation operators for one- 
dimensional Schrodinger operators with potentials, which are asymptotically 
close to Bohr almost periodic infinite-gap potentials. 



1. Introduction 

In the investigation of direct and inverse scattering problems transformation 
operators play an important role. They first showed up in the context of general- 
ized shift operators in the work of Delsarte [5] . They have also been investigated 
by Levitan [25] and were constructed for arbitrary Sturm-Liouville equations by 
Povzner |35j . Afterwards transformation operators have been applied for the first 
time when considering inverse spectral problems, for example by Marchenko |31) . 
who noticed that the spectral function of a Sturm-Liouville operator determines 
the operator uniquely in [30]. Soon after that Gel'fand and Levitan [13] found a 
method of recovering Sturm-Liouville equations from its spectral functions, using 
the transformation operator techniques. 

Another important step was the introduction of transformation operators, which 
preserve the asymptotic behavior of solutions at infinity by Levin |24) . Since that 
these transformation operators are the main tool for solving different kinds of scat- 
tering problems, mainly in the case of constant backgrounds. They have been 
studied for periodic infinite-gap backgrounds by Firsova |10j . and they have 
been recently investigated in the finite-gap case by Boutet de Monvel, Egorova, 
and Teschl 0. 

In the present work we propose a complete investigation of the transformation 
operators for Bohr almost periodic infinite-gap backgrounds, which belong to the 
so-called Levitan class. It should be noted, that this class includes as a special case 
the set of smooth, periodic infinite-gap operators. 
To set the stage, we need 

Hypothesis H. 1.1. Let 

< E a < Ex < ■ ■ ■ < E n < . . . 

be an increasing sequence of points on the real axis which satisfies the following 
conditions: 
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(i) for a certain I > 1, Y^=i{Ein-i) 1 {E^n - E^-i) < oo and 
(ii) E 2n +\ — -E^n-i > Cn a , where C and a are some fixed, positive constants. 

We will call, in what follows, the intervals (2£y_i, E 2 j) for j = 1,2, . . . gaps. In 

each closed gap [£27-1, E 2 j] , j = 1, 2, . . . , we choose a point /x,- and an arbitrary 

sign o-j € {ill- 
Consider next the system of differential equations for the functions fJ-j(x), Cj(x), 

j = 1,2,..., which is an infinite analogue of the well-known Dubrovin equations, 

given by 

(1.1) = - 2a j (x) y /-( H {x) - E )y/ H (x) E 2] ^ H {x) E 2] 

x tt ~ E 2 k-ly/tij(x) - E 2 k 

with initial conditions fJ,j(0) — fij and Oj(0) = o-j, j = 1,2,... 0- Levitan [26], 
[27] . and [55], proved, that this system of differential equations is uniquely solvable, 
that the solutions [lAx), j — 1,2,... are continuously differentiable and satisfy 
fij(x) <G [E 2 j_i, E 2 j] for all x E K. Moreover, these functions Hj(x), j — 1,2,... 
are Bohr almost periodic 0. Using the trace formula (see for example [28] ) 

00 

(1.2) p(x) =E + + E 2] - 2p i {x)), 

j'=i 

we see that also p{x) is real Bohr almost periodic. The operator 

d 2 

(1.3) L := -— +p{x), dom(L) = H 2 (R), 

in L 2 (R), is then called an almost periodic infinite- gap Schrodinger operator of the 
Levitan class. The spectrum of L is purely absolutely continuous and of the form 

a = [Eq, £1] U • • • U [E 2] ,E 2j+1 \ U . . . , 

and has spectral properties analogous to the quasi-periodic finite-gap Schrodinger 
operator. In particular, it is completely defined by the series aj), which 

we call the Dirichlet divisor. Analogously to the finite-gap case this divisor is 
connected to a Riemann surface of infinite genus, which is associated to the function 
Y 1 / 2 {z), where 



(1.4) Y(z) = -(z-E )U 



(z - E 2J ^) (z - E 2J ) 



j=1 E 2 j-i E 2 j-i 

and where the branch cuts are taken along the spectrum. It is known, that the 
Schrodinger equation 

d 2 \ 

— +p(x)Jy(x) = zy(x) 

with any continuous, bounded potential p(x) has two Weyl solutions ip±(z,x) nor- 
malized by 

V>±(z,0) = l and ip±(z t .) E L 2 (R±), for z € C\a. 



^We will use the standard branch cut of the square root in the domain C \ R+ with Im yfz > 0. 
2 For informations about almost periodic functions we refer to |29] , 
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In our case of Bohr almost periodic potentials of the Levitan class, these solu- 
tions have complementary properties similar to properties of the Baker- Akhiczcr 
functions in the finite-gap case. We will briefly discuss them in the next section. 

The objects of interest, for us, are the Jost solutions of the one-dimensional 
Schrodinger operator L in L 2 (R) 

d 2 

(1.5) L:=-—+q(x), dom(i) = H 2 (M), 

with the real-valued potential q(x) £ C(R) satisfying the following condition 

(1.6) / (1 + \x\ 2 )\q(x) ~p(x)\dx < oo. 
Jw 

We will prove the following result. 

Theorem 1.2. Assume Hypothesis Let p(x), defined as in (11.21) . belong to 
the Levitan class, and q(x) satisfy (|1.6[) . then the Jost solutions (f>±(z,x) can be 
represented in the following form 

/*±oo 

<f>±(z,x) = ip±(z,x) ± / K±(x,y)ip±(z,y)dy, 

J X 

where the solutions of (I3.21[) . are real valued, continuously differentiable with respect 
to both parameters and for ±y > ±x they satisfy 

/■±oo 

\K±(x,y)\<±C±(x) \q(s) - p(s)\ds. 

I x+y 

Here C± (x) are continuous positive functions, which are monotonically decreasing 
as x — > ±oo. 

Here it should be pointed out that a subset of the operators belonging to the 
Levitan class consists of operators with periodic potentials. Assume that the se- 
quence {Ej}j° =1 fulfills Hvpothesis ll.il then there is a criterion when this sequence 
is the set of band edges of the spectrum of some Schrodinger operator with periodic 
potential p(x + a) = p(x) > withp £ W% [0, a]. Namely, Marchenko and Ostrovskii 
proved in [32], that 

oo 

P £W*{0,a] iff ^j 2k+2 (y/E^ - Eo - y[E^~ x - E Q ) 2 < oo, 
i=i 

for k = 0, 1, As it is well known that in the periodic case i?2j-i = j 2 + 0(1) 

and E2j — j 2 + 0(1) as j —> oo, we obtain for large j that 

E 2j -i(E 2j - Eij-i) < E 2j _ 1 (^jE 2j -E + y/Eij-i - E )x 

(V E 2 j — E — yj E 2 j-i - Eo) 

< 2j 3 (y/E 2j -E - yfE 2 ~ 1 - E ) =: 

Since Hypothesis 11.11 is satisfied, we have J2jLi J 2k 4 Ij < oo for fc > 2 and hence 
the Cauchy inequality implies that X^Li Ij < oo in this case. This means, that 
Hypothesis 11.11 is satisfied for any a-periodic potential p(x) > with p £ W^fOjo] 
for k > 2. 

The Jost solutions <fi±(z,x) are bounded for z £ a(L) and in particular there is 
no subordinate solution for z in the interior of &(L). Hence the essential spectrum 
is purely absolutely continuous (cf. [15], pi], [H],[l9], [20], and 00] Thm. 9.31]): 



4 



K. GRUNERT 



Corollary 1.3. The essential spectrum of L is purely absolutely continuous 
(1.7) a ac (L)=a{L), a sc {L)=%. 



The point spectrum is confined to the closed gaps M\a(L). 

Criteria that there are only finitely many in each closed gap can be found in [22] , 
[23], US Thm. 6.12]. 

Another application is given by the scattering problem for the operators (|1.3[) 
and (|1.5[) , which has not been considered so far. Such scattering problems arise in 
various physical applications, for example, when studying properties of the alloy 
of two different semi-infinite one-dimensional crystals. A more detailed discussion 
of the history of such problems and further references to the literature can be 
found in [16]. Rather recently the scattering problem has been solved in the case 
of Schrodinger operators with steplike finite-gap potentials in [2] and these results 
have been applied to solve the Korteweg-de Vries equation in [B] and [7J (see also 
[5J GO [33J). Even in the case of periodic backgrounds scattering theory has been 
developed [11] and the time evolution of the scattering data for the Korteweg-de 
Vries equation has been computed [H] . All these investigations use transformation 
operators as a main tool. 

2. Background Schrodinger operators 

In this section we want to summarize some facts for the background Schrodinger 
operator of Levitan class. We present these results, obtained in [55], [37], and [35], 
in a form similar to the finite-gap case used in [2] and [14] . 

Let L be as in (|1.3|) . Denote by s(z, x), c(z, x) the sine- and cosine-type solutions 
of the corresponding equation 

(2.1) (--jjL + p(x)j y(x) = zy{x), z e C, 

associated with the initial conditions 

s(z, 0) = c'(z, 0) = 0, c(z, 0) = s'(z, 0) = 1, 

where prime denotes the derivative with respect to x. Then c(z, x), c'(z, x), s(z, x), 
and s'(z,x) are entire with respect to z. They can be represented in the following 
form _ 

. . , /- s f X siri(\/z(x — y)) . . . 

c{z, x) = cos(^zx) + / -= p(y)c(z,y)dy, 

Jo V z 

, . sm{\fzx') f x sm(^/z(x — y)) . . . . , 
s(z,x) = r ' + VV V y,, p(y) S (z,y)dy. 

v z Jo \ z 
The background Weyl solutions are given by 

(2.2) ip±{z, x) = c(z, x) + m±{z, 0)s(z, x), 
where 

,,,, / x V±(z,x) H(z,x)±Y^(z) 
( 2 -3) m±{z, x) = — = — , 

ip±(z,x) G(z,x) 

are the Weyl functions of L (cf. 28 J, where Y(z) is defined by (11.41) . 
(2.4) G(z : x)^f{ Z ~ E ^ {x \ and H(z,x) = ~G{z,x). 

2 j 1 
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Using (fl~T|) and §2A$ . we have 

The Weyl functions m±(z,a;) are Bohr almost periodic as the following argument 
shows: For each j £ N the functions fJ-j(x) are almost periodic and hence, as a finite 
product of almost periodic functions is again almost periodic, also 



(2.6) G n (z,x) = l[ 



z — p,j (x) 

3 = 1 E V~ l 

is almost periodic for fixed z € C. Moreover, we have 



(2.7) \G n (z,x)-G(z,x)\ 



n / \ oo 

n Z - Hj(X) / _ -r-r 

3=1 j—n+1 J 



and for every fixed zeC there exists a to € N such that \z\ < E^m-i- Then for 
n > to, we obtain on the one hand 

j—n+l J j—n+1 J j=n+l J 

and on the other hand 

n\ Z-Mi x ) \ < TT + N 

j—n+1 J j'=n+l 



(2.9) 



1 00 X 1 

<exp(- Y, (^2i-^--i) + W E ^r- 



where we used that log(l + x) < x and log(l — x) > for x > 0. Noticing 
that the second condition in Hypothesis 11.11 implies that J^JLi x converges, all 
terms are well defined, and it follows that the product XX 00 , , 2 ~^ J ^ converges 

7 * -lxj— n-t-± hj%j_\ ° 

to 1 as n — > oo. Furthermore, Jljli i * s uniformly bounded with respect 
to a; for any fixed z. As all our estimates are independent of x, we have that 
G n (z,x) converges uniformly for fixed z against G(z,x) and thus, the function 
G(z, x) is almost periodic with respect to x. Furthermore by definition g7~~j = 

\ 'g(z 'x) = kO°&(G(z, x)))' and therefore is also almost periodic, where we 

use that log(G(z, x)) ^ for z £ [Ei : E 2 ] U • ■ • U [E 2 j-i,E 2 j] U . . . together with 
[291 Property 3,4,5], and hence m±(z,x) are also almost periodic functions. 

Lemma 2.1. The background Weyl solutions, for ztC, can be represented in the 
following form 

(2-10) ^ 

i, ± {z,x) = exp (I m ± (z,y)dy) = exp (± ^ • 
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If for some e > 0, \z — ^j(x)\ > e for all j € N and x <E M., then the following holds: 
For any 1 > 6 > there exists an R > such that 

(2.11) \tp±(z,x)\ < e T(1 - s)xlm{ ^ ) (l + for any \z\ > R, ±x>0. 

where Dr denotes some constant dependent on R. 
Proof. First we will show that 

G(z,x)^ 1/2 



(2.12) 
fulfills 
(2.13) 



f±(z,x) 



G(z,0) 
d 2 



exp ± 



yV2( z ) 



dx 2 



o G{z,y) 
+ p(x) ) y{x) = zy(x). 



dy 



Using we obtain that f'± = m±/± and = (m' ± +m|)/ ± . Hence ([2~T5|) will 
be satisfied if and only if 

(2.14) m' ± + m± = p(x) - z. 

This can be checked using the following relations, which are proved in [25] , 

(2.15) G(z, x)N(z, x) + H{z, x) 2 = Y(z), 
where 

Z-Tj(x) 



(2.16) 



N(z,x) = -(z-t q (x))YI 



3=1 



Eij -i 



with ro(x) € (— oo, So] and Tj(x) E [J5aj— l, E23] and 

(2.17) ^ G(Z ' x) = 2((p(a;) ~ Z)G(Z ' x) ~ ^ (Z ' 

Moreover, for z outside an s neighborhood of the gaps we have f(z, 0) = 1, and we 
conclude that 



(2.18) 



G(z, x) 



G(z,0) 



n 

3=1 



z — fij [x) 



cxp( £ log(l + *faM) 



3=1 



Thus we obtain 
(2.19) 

G(z, x) 



G(z,0) 



^ ex P ( X! lo S ( 

3=1 



1 



z-fij{0) 



< 



exp 



00 

(E 

3=1 



where we used that log(l + x) < x for x > 0. Moreover 



(2.20) 



Mj(0) ~Mj(g) 



z-/ij(0) 



1 



which implies for |z| < 2i?2j that 
Oi(0) 



(2.21) 



z-Mj(O) 



< 



2( Mj (0)-/i i (x))£ 2j 



THE TRANSFORMATION OPERATOR FOR SCHRODINGER OPERATORS 



7 



For \z\ > 2i?2j we can estimate the terms by 
(2.22) 



( H {0)-^{x))z 



A*i(0) 



< \ N (Q)- N {x)\ 

< \^(0)-^(x)\ 



1 



Z 

1 



< 1^(0) - tij(x)\- 



1 - 



1 - 


E 2j 




z 



<2|/ii(0)-Mi(x)|. 



Combining the estimates from above, we obtain 



(2.23) 



(/ijiO) - fij{x))z 



Mi(0) 



< 2max(l, -^)|^(0) - (j,j(x)\ 

< 2ma X (l,^)(£ 2j - 



and for any fixed e > there exists a fc independent of x and z such that E 
for all n > k, and therefore 



> 1 



(2.24) 



G(z, x) 



G(z,0) 



< exp 



1 



(/ij(0) - Aij(x))z 



Mj (0) 



< exp 



where C is a constant independent of x and z. Analogously one can now investigate 



G(z,x) 

(2.25) 



. Using 



r 1 / 2 (z) = i v ^~^n 

where the roots are defined as follows 
(2.26) 



E 



2j-l 



E- 



2j 



2j-l 



\Tz~~E~ = y/\z - E \e [!iIS< - z - E ^ 2 , 



together with 



y 1/2 i£) : 



is a Herglotz function and 



1 



Eo = Vz(l + O(-)), as z ^ oo, 

z 



(2.27) 

we obtain the following estimate, which is uniform with respect to x, 
Y 1 ' 2 {z) 



(2.28) 



GM =VS(l + 0(-)), asz^oo. 



Using now that J Q X ^Qfj^pj-dr 



G(z,i) 



, where £ € (0, x) by the mean value theo- 



rem, we finally obtain that f±(z, x) has the following asymptotic expansion outside 
a small neighborhood of the gaps asz^oo 



(2.29) 



/±(z,x)=e ±i ^( 1 +°(i))(l + 0(i)), 



where the O(j) terms are uniformly bounded with respect to x and where we use 
the branch cut of the square root in the domain C\R+ with hn(s/z) > 0. Thus 



f±( z , •) <= L 2 (R±) for z € C outside a small neighborhood of the gaps and therefore 
away from the gaps f±(z, x) must coincide with il>±(z, x). This implies in particular 
that f±(z, x) is holomorphic as a function of z in the domain C\{z e C | lm(z) < e}. 
For any point Zo € C with Im(zo) = 2e, we can replace /±(z, x) by its Taylor series 
around the point zo, which has at least radius of convergence e. This Taylor series 



s 
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must coincide with the Taylor expansion of ip±(z,x) around 20, which converges 
on the disc with radius 2e as ip±(z,x) are holomorphic on the domain C\K. Thus 
f±(z, x) is holomorphic on the domain C\R and hence ijj±(z, x) = f±(z, x) on C\K. 
Moreover, this implies that f±(z, x) and ip±(z, x) must coincide on the real axis 
and hence, we finally get that ip±(z,x) can be represented for any z e C by (|2.10l) 
and (|2~29l) implies f2~TT|) . □ 

As the spectrum consists of infinitely many bands, let us cut the complex plane 
along the spectrum a and denote the upper and lower sides of the cuts by a u and a 1 . 
The corresponding points on these cuts will be denoted by z u and z l , respectively. 
In particular, we introduce the notation 

f(z u ) :=lim/(z + ie), f(z l ) := lira/ (2; - ie), z £ a, 
whenever the limits exist. Define the Green function (see e.g. pQ, [4], [21], and [37] ) 

P-30) *,) - -^L, 

where the branch of the square root is chosen in such a way that 

(2.31) tsO") = lm(g{z u )) > for z £ cr, 
1 

then we obtain after a short calculation 

(2.32) W(il>-(z),i/>+(z)) = m+(z) - m-(z) = -g{z)-\ 

where W(f, g)(x) = f(x)g'(x) — f'(x)g(x) denotes the usual Wronskian determi- 
nant. 

For every Dirichlet eigenvalue fij = (J,j(Q), the Weyl functions m±(z) might have 
singularities. If \ij is in the interior of its gap, precisely one Weyl function m + or 
TO_ will have a simple pole. Otherwise, if fj,j sits at an edge, both will have a square 
root singularity. Hence we divide the set of poles accordingly: 

M_|_ = {fij I fij € (E23—1, E2j) and m + has a simple pole}, 

M_ = I Hj G (E23-1, E23) and m_ has a simple pole}, 

M = {Hi I ^ E {E 2j -i,E 23 }}. 

In particular, the following properties of the Weyl solutions are valid (see, e.g. 
[3], [2H], [37], 00]): 

Lemma 2.2. The Weyl solutions have the following properties: 
(i) The functions ip±{z,x) are holomorphic as a function of z in the domain 
C\ (crU M±), real valued on the set R\cr, and have simple poles at the points 
of the set M±. Moreover, they are continuous up to the boundary a u U a 1 
except at the points in M and 



(2.33) V>+0") = ^-{z l ) =^+{z l ), zea. 
For E 6 M the Weyl solutions satisfy 

(2.34) ip±(z,x) = O ( ) , asz^EeM, 

\Vz- EJ 

where the 0((z — E)^ 1 ' 2 ) term is independent of x. 
The same applies to ip'±(z,x). 
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(ii) The functions ip± (z, x) form an orthonormal basis on the spectrum with respect 
to the weight 

(2.35) dp(z) = ^-g(z)dz, 

27T1 

and any f(x) £ L 2 (M.) can be expressed through 

(2.36) f{x) = 1(1 f(y)^+(z,y)dy) ^_(z,x)dp(z). 



Here we use the notation 
(2.37) <f f(z)dp(z) := f f(z)dp(z) - f f(z)d P (z). 

Proof. (i) Having in mind (|2.10l) . we will show as a first step that J* Y q( z [?) dr 

is purely imaginary as z — ¥ E 2 j, with z £ a, (the case z — > E%j-i can be 
handled in the same way). For fixed x £ K we can separate the inter- 
val [0,x] into smaller intervals [0,xi] U [xi,X2] U ■■■ U [xfc,x] such that 
fij(xi) £ {E 2 j-i,E 2 j} and p 3 {x{) ^ p. } {xi +1 ). Assuming fij(xt) = E 2 j-i 
and jj,j{xi+\) = E 2 j 7 and setting 



Y s {z,x) := y/-{z-Eo)^z-E2j-iJ{ 



z - E 2 i-\ z - En x ' 
z - pi(x) z - pi(x) 



where Yj(z,x) is bounded for any z inside the j'th gap, we can conclude 



J xi G{z,t) v J z-^{t) 



dT 



V^e 2 ~( 



z -^ji T ) Jx, z -h(t) 

dT dr+ 



x ' +1 d ~ \ 

Yj(z,T)\ g=Cj ( T )dT), 



dz 



where Cj( T ) G (f J >j( T )> z )- Note that the function -j^Yj(z,T) is uniformly 
bounded for z £ [E^j—i+e, E 2 j+e\ for some e > and that y j( z ^ T ^~^j ^j( r )' r ) 
is uniformly bounded for fij £ [E 2 j—±, E 2 j—\ + e] and z near E 2 j, which 
yields 

(2.38) ^^7/ Tz fj{z,T) ^ {T)dT = ° {^-^ 
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(2.39) 



(2.40) 



(2.41) 



(2.42) 



(2.43) 



On each of the intervals [xj,a;z+i] the function aj(x) is constant and 
therefore 



dr 



E 2j - 



2<Tj - E 2 j (z- fij (t) ) 



-dr 



2j 



= — (Tji arctan 



20]\/y - E 2j {z ~ y) 
1 



dy 



y/E 2j - E 2 j-l \ 



ds 







E 



2j 



A close look shows that the same method can be applied to compute 

JO G(z : t) 



T0j-dr and p x 
This implies that 



k G(z,r) 
xi+i Y 1/2 (z) 
G{z 



dT 



^ydr — >• — \o,jYK as z — > E27 and thus lim^s^ 



G(z,r) 



In more detail one obtains that 



lim exp 

z^E 



x Y^jz) 
o G(z,t) 



dr 



'±1, N {Q)^E, H (x)^E, 
±1, f i j (Q) = E,ii j (x)=E, 
±i, n j (Q) = E, H (x)?E, 
±i, n j (0)^E, f i j {x)=E, 
for any E E da. 

For the product term in (|2.10[) . we have the following estimate 



' exp I ~ a _ ^-i) J - 11 



z-^-(O) 



i=i 



MO) 



z — /ij (x) / 2 

-7^(o) exp v/? 



/2 °° 



for z inside the interval [E 2 j-i + £,E 2 j + e] <~) a for some e > 0, where 
(3 = min Lj . i^j{\E 2j - E 2t \, \E 2 j-i -E 2l \}. This finishes the proof of (j2~Hl) . 
as ([2331) follows directly from ([231]) . 

To show that a similar estimate holds for ip' ± (z,x) as z — > £ 6 M, 
consider, using ([2.101) . 

-0±(z,a;) = m±(z,x)-0±(z,a;)- 

By the investigations of ip±{z, %) from before, it suffices to analyze m±(z, x), 
given by (|2.3[) . which has the following representation 



m±(z,x) = J2~ 



a n (x)Y l / 2 {^ n {x)) 



± 



Y l l 2 (z) 

r[ -f-G(fi n (x),x)(z - /i„(x)) ^ G(z,a:) ' 
For j ^ n we have 



vE^n—Eo 





(2.44) 
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where C\ := exp ^ X^i(-^2j — E2j-i)\ Thus, using Hypothesis 11.11 

we obtain that ^|iMg^ijgL converges uniformly and is uni- 

formly bounded with respect to x. 
For j = n, we obtain 



yV2 (/i . (a;)) 



^G(/i i (a;),x)(z-^(^)) 
and analogously 



. „ y/E 2j - Eoy/in^x) - E 2j _ 1 )(E 2j - Hj(x)) 

— Gl M ■ 



(2.45) 



G{z,x) 



z-Hj(x) 

Multiplying now (|2.43j> . (|2.44|) . and (|2.45p by tp±(z, x), using the estimates 
obtained when proving (|2.34[) , and letting z — > E 2 j with z £ a finishes the 
proof. 

(ii) For a proof we refer to [I Sect. 9.5], [28] Thm 7.1.1], and 41, Sect. 3.8]. 

□ 

3. Derivation of the integral equations for the transformation 

operators and estimates 

Consider the equation 

dx 2 



(3.1) [ -— +q{x)) y(x) = zy{x), z £ C, 



with a potential q(x) satisfying (|1.6|) . Suppose that this equation has two solutions 
4>±{z, x), which we will call the Jost solutions, which, for fixed z, are asymptotically 
close to the background Weyl solutions ip±(z, x) defined in (|2.2|) as x —> ±oo. Set 

/, > T/ \ *P+(z,y)tp-(z,x) - ip + (z,x)ip-(z t y) 

3-2 J (z, x, y) = — -rr—TTT 

W r (V+(«),V'-w) 

and 

(3.3) q(x) = q(x) - p(x). 

It will be shown that the Jost solutions satisfy the following Volterra integral 
equations 

/•±oo 

(3.4) 4>±(z,x) = il;±(z,x)- J(z,x,y)q(y)4>±(z,y)dy. 



We will also show that solutions of this form permit the following representation 

/■±oo 

(3.5) <f>± (z, x) = tp± (z, x) ± / K± (x, y)ij)± (z, y)dy, 

where K±(x,y) are real- valued functions. For simplicity and because both repre- 
sentations can be obtained using the same techniques we will investigate the + case 
only. 

Assume that there exist K+{x,y) with K+{x, .) £ i 2 (R) and K+(x,y) — for 
y < x, such that (f>+(z,x) can be represented by (|3.5|) . Then substituting (|3.5j) 
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into (|3.4|) . multiplying it with ip-(z,x)g(z), integrating over the set a 11,1 , using the 
identity (|2.36j) . and taking into account that K+(x, y) = 0, x > y, we obtain 



(3.6) K + (x,s) + dyq(y) <h J(\,x,y)ilj + (\,y)ip-(\,s)dp(X) 

J X J (7 

poo pOO P 

+ / dyq(y) / dtK+(y,t) f J(X,x, y)tp + (X, t)t/j-(X, s)dp(X) = 0. 

Set 

(3.7) T+(x, y, t,s) = j ^+(A, a:)V-(A, y)4>+(\ t)V-(A, s)g(X)dp(X), 

where the integral has to be understood in the principal value sense 0. 
Then substituting (pL32]l . ([2T35]) . (|3l2]l . and (JBTTJ) into ([H]) we obtain 



(3.8) K + (x,s)+ (T + (x,y,y,s)-T + (y,x,y,s))q(y)dy 

J X 

pOO pOO 

+ d yq(y) / (r+(x,j/,t, s) -r + (y,x,t, s))d* = o. 

A simple calculation using (|2.31l) and (I2.33|) shows that (|3.7[) satishes 



(3.9) T+(a;, y, t, s) = -T+(y, x, s, t). 

Combining (I2.10j) . (I2.30[) . and (|2.35p . one obtains that the only poles of the inte- 
grand are given at the band edges. Using a series of contour integrals and applying 
Jordan's lemma (see for example [39]), it turns out to be necessary to investigate 
the following series 



(3.10) D+(x,y,r, s) = -j f + (E,x,y,r,s), 

where, for any E € da, we denote 



4 

Eeda 



G(z 0) 2 

(3-11) f+(E,x,y,r,s) = lim ' ip+(z,x)il;-(z,y)ijj + (z,r)ijj_(z,s). 

z ^e d_Y(z) 

Lemma 3.1. The series D + (x,y,r, s) defined by A3.10\) and A3.11\) converges, is 
continuous, and uniformly bounded with respect to all variables (x,y,r,s) £ R 4 . 

Proof. First note that by (|2.10[) we have 

(G(E,x)G(E,y)G(E,r)G(E,s)) 1 / 2 



f+(E,x,y,r, s) 



x hm exp / — ' dr + / — > dr 

\J G{z,t) J s G(z,t) 



For informations about principal value integrals, we refer to 1341 Chap. 2] 
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for some E 6 da, where the limit is taken from inside the spectrum. For computing 
the integral terms we refer to the proof of Lemma 12.21 We will now investigate 

(G{E 2l , x)G(E 2l , y)G{E 2l , r)G(E 2l , s)) 1 ' 2 



M + (E 2l ,x,y,r,s) 
(3.12) 



£Y(E 2l ) 

{{E 2l - in{x)){E 2l - My))(E 2 i - Vi{r)){E 2l - M*))) 1/2 



x 



n 



[E 2 i — E )(E 2 i — E 2 i-i) 
((E 2l - iij(x))(E 2 i - ^(y))(E 2l - H (r)){E 2l - Mj (s))) 1/2 



._, (E 2l - E 2 ,j-i){E 2 i - E 2j ) 

which can be estimated as follows 

TT | {E 2 i - Hj(x)) {E 2l - njjy)) < ^ {E 2l - ^(x)) yr {E 2 i - njjy)) 

j= [i# {E21 - E2j) {E2L ~ E *- l] ~ M {E21 ~ E ^ ] iiii {E2L ~ E2 ^ 



.ex P fglog(l + ^-i|^) + £log(l 

\j=i v j=i+i 



t*j(y) -E- 



E 2 j — E 2 i 



1 



DO, 



where we used log(l+x) < x for x > and /3 = min/j^llE^ — \E 2 j-\ — E 2 i\}. 
Moreover, 

(3.13) 

(C^j - w(£)Xg2j - My))(E 2 i - yi{r))(E 2l - Ms))) 1/2 < (E 2l -E2i-!) 

(E 2 i — E )(E 2 i — E 2 i_i) (E 2 i — E ) 

This implies 

(3.14) < \f+(E k ,x,y,r,s)\ = \M+ (E k ,x,y,r,s)\ < d ^ ~ ^ , 

where € {2/ — 1,2/} and Ci := exp f^X^iC^j — E 2 j_i)\ and therefore our 
series converges and hence D + (x, y, r, s) is well defined and uniformly bounded with 
respect to all variables. 

The continuity in all variables follows immediately by using (|2.10[) and (|2.40|) . □ 

Lemma 3.2. The function D + (x,y,r, s), defined through (|3 . 1 0|) and (|3 . 1 1 [) . has 

first partial derivatives, which are uniformly bounded on R 4 . 

Proof. Consider the integral representation (12.101) of the background Weyl solu- 
tions, then the derivative with respect to x is given by 

Note that G(z, 0) 1/2 ^ip+{z, x), by Lemma 12.21 has neither poles nor square root 
singularities at the band edges E G da, which allows us to pass to the limit in the 
following expression 

(3.16) lim (z - E) H[Z ' X) + yl/2(z) (G(Z ' y)G{z > r)G(Zj S))1/2 



z->E y ' G 1 / 2 {z,x) Y(z) 
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Therefore, we will slightly abuse notation by omitting the limit and replacing z by 
E. W.l.o.g. we will assume that E = E2n (the case E = E 2n -i can be treated 
similarly). Using (12.51) . we have 

H(E 2n ,x) (G{E 2n ,y)G(E 2n ,r)G(E 2n , S )) 1 / 2 



GV*(E 2n ,x) 



£Y(E 2n ) 



M+(E,x,y,r, s). 



Due to (13.141) we will at first consider V\ , d Y ^jf^ — rrr- 
done using the same ideas as in Lemma 13. II Namely, for j =/= n, 



which can be 



(3.17) 0< 



YV*{ H {x)) 



4-G(nj{x),x)(E 2n -Hj{x)) 



1 



< -p\l{E 2j - Eq){E 2j - E 2j -i)Cl /2 : 



which implies that the corresponding sum converges. 
For j = n, there are two cases to distinguish: 



(i) If (i n (x) = E 2n , then 
\Y^(E 2n " 



< 



G(E 2n ,x) 



M + (E 2n ,x,y,r,s) 



3/2 ((E 2n - Hn(y))(E 2n ~ Hn(r))(E 2n - ^„(s))) 1/2 



\J{E 2n — E )(E 2n — E 2n -\) 



< £,3/2 E 2n — E 2 n-1 

V E 2n — Eq 
(ii) If fi n (x) ^ E 2n , we have 

YV*{» n {x)) 



< 



-j^G(n n (x),x)(E 2n - fi n (x)) 



M + (E 2n ,x,y,r, s) 



n 3 / 2 f I \ TP M/2 ((Mg) ~ E 2n-l){E 2 n - /i„(?/))) 1/2 

<C 1 {fJ.„{X) - E ) i — — 

\&2n — ) 

{{E 2n - n n (r))(E 2n - Mn(s))) 1 / 2 



(E 2n — E 2n -i) 



Next we consider 



3/2 (E 2n — E 2n -i) 



(3.18) 



V E 2n — Eq 



Y Y l 2 {E 2n ) (G(E 2n ,y)G(E 2n , r)G(E 2n , s)) 1 / 2 



GV 2 (E 2n ,x) 
which can be investigated as before. 



1-Y(E 



2n ) 
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(i) If fi n (x) = E 2n , then 



< 



G(E 2n ,x) 



M+(E 2n ,x,y,r,s) 



< c ,3/2 (Cggn - Vn(y))(E 2n - Vn(r))(E 2n - fi n (s))) 1/2 

\J (E 2n — E )(E 2n — _E 2 n-i) 

^ r ,3/2 {E 2n — E 2n -l) 
- °1 



\J (E 2n — Eo) 
(ii) If li„(x) ^ E 2n , 

Y^(E 2n ) 



G(E 2n ,x) 

because 



M + (E 2n ,x,y,r, s) = 0, 



/ 3 ig) Hm (-( z - E 2n -i)(z - E 2n )(z - fin(y))(z - fin(r))(z - ^n(a))) 1/2 _ Q 

«-^^» y/(z - E )(z - i*n(x))(z - £ 2 „-i) 

This finishes the proof, as all our estimates are uniformly with respect to x, y,r, s £ 
K and 

(3-20) 0<|/+,,(£ fe , WlS )|<^%=i%^, 

V (^2n — J^oJ 

where k £ {2n — 1, 2n} and C denotes a constant independent of n. □ 

Lemma 3.3. The kernels K±(x,s) of the transformation operators satisfy the in- 
tegral equation 

/■±oo 

K±(x,s) = -2 / q(y)D±(x,y,y,s)dy 

/*ioo ps-\-y—x 

(3-21) T2/ dW I>±(a!,y,r,8)lir ± (»,r)g(y)dr, ±s > ±s, 

where D±(x,y,r, s) are defined by h3.10\) . In particular, 

(3.22) K ± {x,x)=±^J (q(s)- P (s))ds. 

Proof. Suppose that (x — y + r — s) > 0, where x, y, r, s are considered fixed pa- 
rameters, and take a series of closed contours C n consisting of a circular arc R n 
centered at the origin with radius (E 2n + E 2n+ i)/2, together with some parts wrap- 
ping around each of the first n + 1 bands of the spectrum a, but not intersecting 
it, as indicated in figure [TJ Then we can consider the following series of contour 
integrals 

(3.23) I n (x,y,r,s) := — [ %p + {z 1 x)^^(z 1 y)i} + (z 1 r)il;^(z,s)g(z)' 2 dz. 

27ri Jc n 

Moreover, consider L n the integral along the loop around the n'th band of the 
spectrum and set 

1 



(3.24) S n = —(b if) + (z,x)<if)_(z,y)i/} + (z,r)il;-(z,s)g(zydz, 

ZlT1 J[E2n-2,E 2 n--L] 
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Figure 1. Contour C, 



then we have 

(3.25) L n = L n S n + S n = S n , 

where we used that L n — S n = 0, by the continuity of the integrand up to the real 
axis and Cauchy's theorem. 

Thus, shrinking the loops around each band of the spectrum in C n to an integral 
along each band of the spectrum in C n , and denoting by I n (x, y, r, s) the resulting 
series of contour integrals, we obtain 

(3.26) I n (x,y,r,s) = / tp+(z, x)ip-(z, y)4>+(z, r)tp-(z, s)g(z) 2 dz 

^ 2n 

- -^^2f+(Ei,x,y,r,s), 

1=0 

where we used Sokhotsky-Plemelj's formula (see e.g. [Ml Chap. 2]). 

On the circle R n we have the following asymptotic behavior as n — > oo and 
therefore z — > oo, 

(3.27) g 1/2 (z)1>+(z, x) = e ^(i+o(i)) ^J_^) , 

where we used Lemma 12.11 and the fact that these asymptotics are valid as long as 
we are outside a small neighborhood of the gaps. This yields 

(3.28) lp+(z, x)tp-(z, y)ip+(z,r)il)-(z,8)g(z) 2 = ^VS(x-y+r-s)(l+0(l)) q 
as z — > oo. 

Hence one can apply Jordan's lemma to conclude that the contribution of the 
circle R n vanishes as n — > oo. 

Thus letting n — > oo the series of integrals I n (x,y,r, s) converges to 

(3.29) T + (x,y,r,s) = D + (x,y,r,s), for (x - y + r - s) > 0. 

Note that f+(E, x, y, r, s) is real for any E £ da, because (|2.40[) and G(E, x) = 
0, if fij(x) = E, imply that f+(E,x,y,r, s) = 0. Moreover f+(E,x,y,r,s) — 
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f+{E, y, x, s, r). Thus D+(x, y, r, s) is also real and satisfies 

(3.30) D + (x,y,r,s) = D + (y,x,s,r). 

Now let {x - y + r - s) < 0, that is -(a; - y + r - s) > 0. Then IpTS]) . (|3~29| . 
and (|3.30[) imply 



(3.31) T + (x,y,t,s) = -T + (y,x,s,r) = -D + (x,y,r,s) = -D + (x,y,r, s). 
Therefore, 

(3.32) T + (x, y, r, s) = D + (x, y, r, s) sign(a; -y + r-s). 

Combining all the facts, the domain, where the first integrand in (|3.8p does not 
vanish, is given by 

(3.33) sign(a; — s) — — sign(2y — x — s), s > x. 
In the second integral the domain of integration is 

(3.34) sign(x — y + t — s) = — sign(y — x + t — s), with s > x, t > y > x. 
Solving (pT33l) and (pT34"l) . proves (f3~2T]) . 

Setting now s — x in (|3.21[) . the second summand vanishes, because we set 
K + (y,r) = for r < y. Hence 



/•OO 

(3.35) K+(x,x) = -2 q(y)D+( 

J X 



x,y,y,x)dy. 

J X 

Thus we obtain 

(z - Hj(x)){z - Hi(y)) 



1 1 

(3.36) D+ (x, y, y, x) = - Res — — JJ 



Be3a - z ~ E ofJ i {z - E 2j -i){z - E 2j ) 



4 .f— ^ e 

z-Ei yt{z- (J.j(x))(z - Hj{y)) 



= - V li z ~ — TT 



2j, 



and we already know that this function is bounded by Lemma 13.11 Considering 
now the following sequence 

(3.37) D +A x,y,y,x)= 1 -j: Res * ft ( f ~ ^Pfe % (g » 



I V lim Z ~ jE ' TT ^-/^(^x^-MjfeQ) 

4 i ^ Z - ^0 f=l ( Z - 3W-l)(* - ^) ' 
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which corresponds to the case where we only have n gaps and crossed out all the 

other ones. We now estimate 

(3.38) 



^ 2n 

\D+(x,y,y,x) - D ni+ (x,y,y,x)\ <-V lim 

4 L — J z^tEi 
1=0 



z-Ei 



z - fij(x) z - Hj(y) 



z- E n y z - E 2 j-i z-E 2 j 



n 

tj—n+l 



^ .X. 

— > lim 

4 ^ z ^ El 

l= 2 n+l 



z - fi^x) z - Hj(y) 
z — E%j-i z — E 2 j 



- 1 



z-E 



n 



z - Hj{x) z - nj(y) 



z-E- 



2J-1 



z-E- 



2j 



using the same techniques as in the proof of Lemma 13.11 We fix z = E% (the case 
z = E 2 i-\ can be handled analogously). If I < n, we have 



lim 

r->_B 2i 



z - E% 



Z-En 



z - fij(x) z - fXj(y) 
l J- i \z — E 2 j-i z — E 2 j 



n 



E 2 i - Hi(x) E 2 i - m{y) 



n 



E 2l - fXj(x) E 2 i - fij{y) 



E 2 i — Eq E 2 i — E 2 i-i . , J r E 2 i — E 2 j-\ E 2 i — E 2 j 



i-i 



< 



< 



n 



E 2t - fij(x) 



Eqi — E 2 l-l y-r E 2 i — fij(x) 

E 2 i — Eq f-f- E 2 i — E 2 j .- L ,- L 1 E% — E 2 j-i 

j = l J 3=1+1 ■> 



and 
(3.39) 



1 X 

exp ( - b Z ( E v 



E 2j-1, 



< 



j=n+l 



n 

j=n+l 



(ij(x) - E 2 i 



E 2 i — E- 



2 j - rj 2i 



< 



n 

j=n+l 



fijjx) - E 2i [ij\y] - E 2X 
E 2 j-i — E 2 i E 2 j — E 2 i 



< 



fijjx) - E 2 i 
Eo<i—-\ — E 2 i 

j=n+l ^ 1 M 



n 



^ 6x p(^ E ( E ^~E 2j -i) 



j=n+l 



where the last estimate implies that the first sequence in (|3.38[) converges uniformly 
to zero as n tends to oo as in the investigation of (|2.7[) . Analogously, one can 
estimate the terms of the second sequence to obtain 



lim 



z-Ei 



n 



z - fij(x) z - Hj{y) 



z — Eo \z - E 2 j-i z - E 2j 



E 2 [ — E 2 i-i v-v, p „ \ 

" E 2l -Ec «P(flE(^-^-0 



3 = 1 



Thus also the second sequence in (|3.38p converges uniformly to zero, as we are 
working in the Levitan class and hence D n ^ + (x, y, y, x) converges uniformly against 
D+(x,y,y,x). 
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Moreover, it is known (see e.g. [2]), that 

(3-40) D n ,+ (x,y,y,x) = -- 

for each fixed n, and hence we finally obtain 

(3.41) D+(x,y,y,x) = lim D n , + (x,y,y,x) = lim -\ = -\. 

Therefore we can now conclude, using (|3.35[) . that 



(3.42) 



K±(x,x) = ±- 



±00 



{q(s) -p(s))ds 



□ 



Lemma 3.4. Suppose \1.6\) . then \3.21\) has a unique solution K±(x,y) such that 
K±(x, y) has first-order partial derivatives with respect to both variables. Moreover, 
for ij/ > ±x the following estimates are valid 



(3.43) 
(3.44) 
where 
(3.45) 



\K ± (x,y)\<C±(x)Q ± (x + y), 



dK±(x,y) 



dx 



dK±(x,y) 



dy 



< C±{x) 



x + y 



+ Q±{x + y)\, 



Q±(x)=± \q(s)\ds, q{x) = q(x) -p(x), 



and C±(x) are positive continuous functions for x € which decrease as x — > ±oo 



and depend on the corresponding background data and on ± J 2 ±oc Q±(s)ds. 

Proof. Using the method of successive approximation one can prove existence and 
uniqueness of the solution K±(x 7 y) of (|3.21l) . We restrict our considerations to the 
+ case. After the following change of variables 

(3.46) 2a := s + r, 2/3 := r — s, 2u := x + y, 2v := y — x, 

(|3.21[) becomes 

r°° 

H(u 7 v) = — 2 / q(s)Di(u,v, s)ds 

J u 

q(a - (3)D 2 {u, v, a, P)H{a, (3)d(3, 



(3.47) 



ith 



da 



H(u, v) — K + {u — v, u + v), Di(u, v, s) = D + {u — v, s, s,u + v), 

(3.48) D 2 {u,v,a,{3) = D + (u - v, a - f3, a + f3,u + v). 

As the functions D\ and D 2 are uniformly bounded with respect to all their variables 
by a constant C, we can apply the method of successive approximation to estimate 
H(u,v), which yields 

(3.49) \H(u,v)\<C(u-v)Q + {2u), 
where 



(3.50) C(u- v) = 2Cexp [AC [ Q^ 

V J2u-2v 



(x)dx 
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To obtain the second estimate p.44j) we recall that the partial derivatives with 
respect to all variables exist for D\ and D 2 and that they are also bounded with 
respect to all variables. Thus, using 

dH(u, v) 



(3.51) 



du 



— 2q(u)Di(u, v, u) 



- i I q(u - (3)D 2 (u, v, u, (3)H(u, (3)d(3 - 2 £ q(s) dDl ^ V,s) ds 
da J q(a-p) — H{a,(3)d/3, 



and 
(3.52) 



dH(u, v) 
dv 

= -2 



1 dDi(u,v,s) 
q(s) - ds + 2 



dv 

da / q(a — (3) 



q(a — P)D 2 (u, v, a, v)H{a 1 v)da 



dD 2 (u, v, a, (3) 
dv 



H 



one obtains 
(3.53) 



\—H(u,v)\ < Ci(u-«)(|?(u)| +Q+(2u)), 

\^-H(u,v)\ < d(u - v)(\q(u)\ + Q + (2u)), 
dv 



where C\(u — v) is of the same type as (|3.50|) with a different constant C\ depending 
on the background data. Using 



(3.54) 



dK+(x, y) 



dx 



dK + (x,y) 



All 



< 



dH(u, v) 



du 



dH(u, v) 



dv 



completes the proof. 

To finish the proof of Theorem ll.21 we have to show the following: 
Lemma 3.5. The functions 



□ 



(3.55) 



4>±{z,x) = ip±(z,x) ± 



±00 



K±(x, s)ip± (z, s)ds, 



where K±{x,s) are defined by (|3.21[) . satisfy 
(3.56) ' ''' 



( - j~2 + q( x ))<t>±(z,x) = z(j)(z,x). 



Proof. Again we will only consider the + case as the other one can be treated 
similarly and drop the + whenever possible. On the one hand we obtain, using 
(|332|) . that 

( d 2 \ 

\ ~ ^2 + q ( x ^) 0( z > x ) = ~^ ( z > x ) + Pi x )^i z i x ) 



(3.57) 



+ -q(x)ip(z, x) + K(x, x)?p r (z, x) + K x (x, x)ip(z, x) 
+ / (q(x)K(x,s) - K xx (x,s))tfj(z,s)ds, 
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and on the other hand, using that ip±(z, x) are the background Weyl solutions, we 
have 

Z(j)(z, x) = zip(z, x) + K(x, x)ip'(z, x) — K y (x, x)ip(z, x) 

(3.58) + / (p(s)K(x,s) - K ss (x,s))ip(z,s)ds. 

J X 

Applying f)3.22jl once more, we see that (|3.56p is satisfied if and only if 

poo poo 

(3.59) / (K xx (x,s) - K ss (x,s))tp(z,s)ds = (q(x) - p(s))K(x,s)ip(z,s)ds. 



For proving this identity we use the integral equation (13.6[) instead of (|3.2ip for 
K{x, s), which yields for x < s 



K ss (x,s) = - / dyq(y) i> J(X,x,y)tp + (X,y)ip"(X,s)dp(X) 

x J a 

I poo p 

dyq{y) / dtK(y, t) f J(A, x, y)i/> + (A, €)^"_ (A, s)dp(X) 

J y J a 

(p(s)-p(x))K(x,s) - / dyq(y) d> J XX (X, x, y)tp + (A, y)tp- (A, s)dp(X) 

J x J a 

i pOO p 

dyq(y) / dtK(y,t) <b J xx {X,x,y)<i/; + (X,i)ip-(X,s)dp{X) 

J y J a 

= (jp{s) -p(x))K(x,s) + K XjX (x,s) - q(x) <j> J x (X,x,x)?p + (X,x)ip-(X,s)dp(X) 
— q(x) / dtK(x,t) (f) J x (X,x,x)ifj + (X,t)il;^(X,s)dp(X) 

J X J (J 

(3.60) 

= (p{s) - q(x))K(x, s) + K xx (x, s). 

Here it should be noticed that <f J(X, x, y)ijj+ (A, t)ip'_!_ (A, s)dp(X) exists, because we 
can again estimate the sum of the absolute values of the residues by using that 
—ip±(z,x) + p(x)ip±(z, x) — zip±(z,x) and the same techniques as in Lemma |3. II 
Analogously for § J XX {X, x, y)tp + (X 7 t)?/>-(A, s)dp(X). Thus (|3.59[) is fulfilled and 
therefore also (|3.56|) . 

□ 

Lemma 3.6. One has K±(x,y) 6 L 2 (M.±,dy). 
Proof. Using (|1.6[) . we conclude 

(3.61) \K ± (x,y)\ 2 dy<±C ± (x) 2 Q±(x + y) 2 dy 

JU Jx 

/•±oo /-ioo 

= C±(x) 2 Q±(2x) / / \q(s)\dsdy 

Jx 

/>±00 

= C±{x) 2 Q ± (2x) / (2s- 2x)\q(s)\ds < oo. 

J X 

□ 
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It should also be noticed that for any function f±(x) € i 2 (M±), 

/■±oo 

(3.62) h±(x) = f±(x)± K ± (x,y)f ± (y)dyeL 2 (R ± ). 

J X 

Thus as a consequence we obtain 

Corollary 3.7. The normalized Jost solutions 4>±{z,x) coincide with the Weyl 
solutions of the Schrodinger operator hl.5\) . 

Acknowledgments. I want to thank Iryna Egorova and Gerald Teschl for many 
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